The probability distribution of the conductance at the mobility edge 
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The probability distribution of the conductance p(g) of disordered 2d and 3d systems 
is calculated by transfer matrix techniques. As expected, p(g) is Gaussian for extended 
states while for localized states it is log-normal. We find that at the mobility edge p(g) 
is highly asymmetric and universal. 
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In the presence of disorderly, a system may undergo 
a transition from insulating to metallic behaviour as 
the Fermi energy varies in an energy range containing 
both localized and extended states, seperated by a mo- 
bility edge. This transition can be characterized cither 
by transport properties like e.g. the conductance, or by 
properties of the eigenstates of the system like e.g. the 
correlation length £ c (approaching from the metallic side 
of the transition) or the localization length (approach- 
ing from the insulating side of the transition). While 
the latter are self-averaging quantities, i.e. the ensemble 
average majhbe use d as a scaling variable, the conduc- 
tance is notnQ. Therefore it is of great importance to 
determine the complete probability distribution p(g) of 
the conductance g (in units of e 2 /h), especially at the 
critical point of the metal-insulator transition, as it is 
well known to be a Gaussian on the metallic side and 
log-normal on the insulating side. The correct form of 
Pc(gk jthe distribution of the conductance at the mobility 
edgeaQ, is still not sufficiently well-known. 

Using a tight-binding modelcl with diagonal disorder, 
a transition from a metallic state to an insulating one 
can be inducedEI in a finite size sample by increasing the 
disorder strength W. In all our results, W is given in 
units of the hopping integral. The localization length 
decreases as the strength of the disorder, W, increases. 
As long as is much bigger than the system size, the 
electron will cross the sample with ease, thus being essen- 
tially delocalized. If, on the other hand, £j is sufficiently 
smaller than the system size, the electron will become 
localized in a small region and not contribute much to 
the conductance. The critical strength of disorder W c 
will occur where the localization length becomes compa- 
rable to the system size. We have systematically studied 
the conductance g of the 2d and 3d tight-binding model 
by using the transfer matrix technique, which relates g 
with the transmission matrix t by g = 2Tr(t T t). The g 
defined here is for both spin orientations. Fig. || shows 
the distributions of the conductance g and the natural 
logarithm of the conductance ln(g) for a metallic and a 
localized 2d sample respectively. In the "metallic" regime 
(it is really weakly localized, since we are working in 2d) 
where W = 3.0 we have £; = 47234, which is much larger 



than the system size L = 64. In the localized regime, 
where W = 9.0, we have = 7.54, which is smaller 
than the size of the system. Both distributions gan be 
fitted very well by a Gaussian normal distribution!! Fig. 
[I] displays the distributions for a system near the critical 
point. The distribution of ln(g) shows a characteristic 
cut-off at ln(g) = 0. Only a few samples have a con- 
ductance g > 1. This is in agreement with analytical 
results!. 




FIG. 1. The distribution of the conductance for a square 
system of 64 by 64 lattice sites. The disorder was chosen such 
that the localization length is close to the system size L. 
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very similar to the case shown in Fig. |]. 
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FIG. 2. The distribution of the conductance for a square 
of 64 by 64 lattice sites. Small disorder (W = 3, top panel) 
leads to metallic behaviour; strong disorder (W = 9, bottom 
panel) leads to insulating behaviour. The values for the aver- 
age and standard deviation of the Gaussian fits are given in 
the graphs. 

Applying a strong magnetic field perpendicular to a 
2d system creates states with a diverging localization 
lengthB even in the thermodynamical limit, as long as 
the disorder is not too strong. Thus, one can approach a 
critical state in such a system by varying the energy, even 
though one cannot reach a true metallic state. In Fig. || 
we have an insulating system, the conductance distribu- 
tion again fitting well to a log-normal distribution. The 
distribution for the critical state p c {g), shown in Fig. ||, 
again has the abrupt cut-off at ln(g) = 0. It can be fitted 
to a skewed log-normal distribution, which is normalized 
on the interval (— oo; 0] rather than all real numbers. No- 
tice that in this case too p c (g) is highly asymmetric and 
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FIG. 3. The conductance distribution of a square system 
of 64 by 64 lattice sites. For the chosen disorder strength 
W = 4 and magnetic flux a = 1/8 states at the energy 
E — —3.8 are well localized. The parameters for the Gaussian 
fit are given in the graph. 



0.8 







2.0 






1.5 




0.6 


"S i.o 








/ — s 




0.5 


So 






c 


0.4 


- 0.0 


ft, 




C 




0.2 





<\n(g)> 



0.0 




ln(g) 



FIG. 4. The conductance distribution of a square system 
of 64 by 64 lattice sites. For the chosen disorder strength 
W = 4 and magnetic flux a — 1/8 states at the energy 
E — —3.4 are critical. The parameters for the fit are given in 
the graph. 

Comparing these results to a 3d system without mag- 
netic field, again varying the disorder strength and keep- 
ing the energy fixed at E = 0.0, we find the same qual- 
itative picture: the distribution of the conductance is 
normal on the metallic side and log-normal on the in- 
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sulating side of the transition (see Fig. ||), whereas the 
critical state is characterized by a cut-off at \n(g) = 
and a skewed log-normal distribution for \n(g) < (see 
Fig. |). 
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FIG. 6. The conductance distribution of a cubic system 
of 10 by 10 by 10 lattice sites. The disorder was chosen such 
that the localization length is close to the system size L. The 
parameters for the fit are given in the graph. 
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In conclusion, our detailed numerical results show that 
the probability distribution of the conductance is normal 
for the extended regime and log-normal for the localized 
regime. However, at the mobility edge p(g) is highly 
asymmetric. The form of p c (g) at the critical point is 
independent of the dimensionality of the system and of 
the model. This suggests that p c (g) is universal. 

Ames Laboratory is operated for the U.S. Department 
of Energy by Iowa State University under Contract No. 
W-7405-Eng-82. 
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FIG. 5. The conductance distribution of a cubic system 
of 10 by 10 by 10 lattice sites. Small disorder (W = 3, left 
panel) leads to metallic behaviour; strong disorder (W = 30, 
right panel) leads to insulating behaviour. The values for the 
average and standard deviation of the Gaussian fits are given 
in the graphs. 
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